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ABSTRACT 

Using  new  algorithms,  we  derive  estimates  of  partitioned  absorption  (a)  and  scattering  (b) 
coefficients  from  high-resolution,  hyperspectral  aircraft  imagery  collected  with  the 
PHILLS  sensor.  During  May  2002,  we  conducted  simultaneous  ship  and  aircraft  surveys 
in  a  dynamic  environment  in  the  northern  Gulf  of  Mexico  (coastal  waters  near  Mobile 
Bay,  Alabama).  Optical  measurements  included  hyperspectral  remote  sensing  reflectance 
(Rrs),  partitioned  absorption  coefficients  (phytoplankton,  detrital,  and  colored  dissolved 
organic  matter  components),  beam  attenuation  (c),  backscattering  (bb)  coefficient,  HPLC 
phytoplankton  pigment  concentrations,  and  organic/inorganic  particle  loads.  Our 
objectives  are  to:  (1)  develop  new  algorithms  to  estimate  optical  properties  from 
hyperspectral  ocean  color  imagery;  (2)  apply  the  algorithms  to  high-resolution  aircraft 
imagery  to  characterize  the  spatial  optical  variability  in  an  environment  impacted  by  high 
concentrations  of  dissolved  and  particulate  materials;  and  (3)  perform  an  optical  water 
mass  classification  on  the  aircraft  imagery.  The  new  optical  classification  system  based 
on  the  partitioned  absorption  coefficients  is  used  to  distinguish  water  masses  and  track 
coastal  features. 

The  PHILLS  radiance  values  are  first  converted  to  atmospherically-corrected  reflectances 
through  an  empirical  line  calibration  technique  that  couples  the  imagery  with  in  situ 
reflectance  measurements.  Then,  by  focusing  initially  on  the  near-infrared  (NIR)  region 
of  the  reflectance  spectrum  where  absorption  is  dominated  by  pure  water,  we  can 
estimate  the  backscattering  coefficient  at  715  nm  by  solving  a  quadratic  equation.  A 
linear  model  extends  this  estimate  spectrally  to  shorter  wavelengths.  From  reflectance 
and  backscattering  we  subsequently  estimate  the  total  absorption  coefficient.  Empirical 
relationships  based  on  filter  pad  measurements  are  used  to  partition  the  total  absorption 
into  phytoplankton,  detrital,  and  CDOM  components,  as  well  as  organic  and  inorganic 
particulate  loads.  Finally,  we  compare  the  aircraft-derived  estimates  with  in  situ  optical 
measurements  collected  during  the  coincident  ship  surveys,  and  assess  errors. 

INTRODUCTION 

Remote  sensing  imagery  can  provide  synoptic  snapshots  of  water  optical  properties. 
Airborne  imagery  generally  provides  increased  spectral  and  spatial  resolution  compared 
to  satellite  imagery  (Hyperion  is  the  only  hyperspectral  satellite  sensor  currently  on  orbit 
(220  channels,  30  m  ground  resolution,  16-day  repeat  cycle),  but  it  was  launched  by 
NASA  as  a  demonstration  project  with  a  one-year  on-orbit  design  lifetime,  and  it  does 
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not  provide  operational  data).  In  addition,  the  temporal  resolution  of  the  aircraft  imagery 
(repeat  coverage)  is  only  limited  by  flight  requirements  and  cost  constraints,  not  by 
orbital  paths  and  times.  Thus,  shorter  temporal  and  spatial  scales  (i.e.,  higher  frequency 
processes  such  as  tidal  variability)  can  generally  be  resolved  with  aircraft  imagery. 

However,  aircraft  imagery  is  also  fraught  with  problems  unique  to  that  platform.  For 
example,  because  the  sensor  is  flown  at  variable  heights  within  the  atmospheric  layer, 
rather  than  above  it  as  is  the  case  for  satellite  sensors,  the  atmospheric  correction  problem 
is  more  difficult,  particularly  over  coastal  areas.  The  platform  itself  is  also  less  stable 
than  a  satellite;  the  aircraft  data  require  roll,  pitch,  and  yaw  corrections  for  accurate 
ground  pixel  locations.  Instrument  calibration  is  an  ongoing  issue  for  any  optical  sensor 
that  is  operated  in  a  harsh  environment,  such  as  the  constantly  vibrating  and  jarring 
conditions  within  an  airplane.  These  difficulties  represent  challenges  that  must  be 
overcome  before  high-quality  data  are  available  for  incorporation  into  bio-optical 
algorithms. 

OBJECTIVES 

Our  objectives  are  to:  (1)  develop  new  algorithms  to  estimate  optical  properties  from 
hyperspectral  ocean  color  imagery;  (2)  apply  the  algorithms  to  high-resolution  aircraft 
imagery  to  characterize  the  spatial  optical  variability  in  an  environment  impacted  by  high 
concentrations  of  dissolved  and  particulate  materials;  and  (3)  perform  an  optical  water 
mass  classification  on  the  aircraft  imagery. 

METHODS 

In  May  2002,  we  conducted  field  experiments  near  Mobile  Bay,  Alabama,  to  characterize 
the  coastal  optical  properties  using  ship,  satellite,  and  aircraft  data.  A  total  of  63  stations 
were  occupied  by  two  ships  over  an  eight-day  period.  One  ship,  the  R/V  Pelican,  mapped 
the  outflow  plume  from  Mobile  Bay  while  the  other  vessel,  the  R/V  Ocean  Color, 
examined  spatial  variability  of  optical  properties  in  Mobile  Bay  and  Mississippi  Sound 
(shoreward  of  the  barrier  islands).  In  situ  optical  measurements  included  hyperspectral 
remote  sensing  reflectance  (ASD  and  Spectrex  instruments),  absorption  and  beam 
attenuation  coefficients  (and  scattering  coefficient  by  difference;  ac9  instrument), 
backscattering  coefficient  (Hydroscat  instrument). 

In  addition,  whole  water  samples  were  collected  on  both  ships.  Filter  pad  absorption 
analyses  were  performed  using  the  ASD  fiber  optic  spectroradiometer  to  partition  the 
total  absorption  coefficient  into  phytoplankton  (a®),  detrital  (a^),  and  colored  dissolved 
organic  matter  (gelbstoff,  ag)  components.  GF/F  filtered  samples  were  also  collected  and 
frozen  in  liquid  nitrogen  for  subsequent  fluorometric  and  HPLC  pigment  analysis.  Total 
suspended  solids  (TSS),  particulate  inorganic  matter  (PIM),  and  particulate  organic 
matter  (POM)  were  determined  gravimetrically.  Prior  to  the  cruise,  47  mm  GF/F  filters 
were  rinsed  three  times  with  Milli-Q  water,  dried  for  two  hours  at  103°  C,  ashed  for  15 
minutes  at  550°  C,  weighed  on  an  analytical  scale,  and  stored  for  analysis  in  small, 
labeled  plastic  petri  dishes.  During  the  experiment,  whole  water  samples  were 
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subsequently  filtered  through  the  prepared  filters,  rinsed  thoroughly,  dried,  and  weighed 
to  determine  TSS  in  mg/1.  The  filters  were  then  ashed  and  re-weighed  to  determine  PIM 
(and  POM  through  subtraction). 

Coincident  aircraft  imagery  was  collected  with  the  Portable  Hyperspectral  Imager  for 
Low-Light  Spectroscopy  (PHILLS).  Four  flight  days  were  successfully  completed  with 
4-8  overlapping  flight  lines  per  day  oriented  parallel  to  the  coast  (into  or  away  from  the 
sun)  to  allow  image  mosaicking  with  minimal  sun  glint.  PHILLS  operated  with  220 
spectral  bands  at  a  flight  altitude  of  10,000  m  to  yield  a  10  m  ground  pixel  spatial 
dimension.  Figure  1  shows  an  example  of  the  flight  lines  and  times,  as  well  as  the 
coincident  ship  stations  and  collection  times,  for  21  May,  2002.  Similar  figures  for 
flights  on  19,  20,  and  23  May  are  not  shown. 
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Figure  1.  Locations  and  times  of  flight  lines  and  ship  stations,  21  May,  2002. 


The  geolocated,  calibrated  radiance  data  (Kohler,  2002)  were  subsequently  converted  to 
remote-sensing  reflectancs  (Rrs)  using  an  empirical  line  calibration  (ELC)  technique 
(Moran  et  al.,  2001).  This  approach  simultaneously  applies  an  atmospheric  correction  to 
the  imagery  by  forcing  the  aircraft  radiance  data  to  match  the  in  situ  Rrs  data.  Following 
the  ELC,  a  near-infrared  (NIR)  slope  algorithm  was  applied  to  the  imagery  to  derive 
estimates  of  inherent  optical  properties  and  chlorophyll.  The  NIR  algorithm  is  described 
below.  Following  application  of  the  NIR  algorithm,  the  imagery  was  goecorrected  and 
mosaicked.  Optical  property  values  were  then  extracted  from  the  imagery  at  the  locations 
of  the  ship  stations  and  compared  to  the  measured  values.  The  steps  are  outlined  in 
Figure  2. 
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Figure  2.  Schematic  of processing  steps  for  PHILLS  aircraft  imagery. 


The  ELC  technique  applies  constant  calibration  coefficients  across  the  flight  line  (and  in 
this  case  across  multiple  flight  lines)  and  thus  may  not  be  as  desirable  as  an  atmospheric 
model-based  correction  scheme  such  as  ATREM  or  TAAFKA  (Gao  et  al.,  2000). 
However,  the  model-based  corrections  require  some  knowledge  (or  estimate)  of  several 
atmospheric  properties,  such  as  ozone  thickness,  aerosol  optical  depth,  relative  humidity, 
etc.,  and  our  experience  with  those  approaches  is  that  they  often  yield  unreliable  results 
and/or  negative  radiances  in  coastal  waters,  so  we  explored  the  ELC  approach. 

As  mentioned  above,  we  used  a  single  station/flight  line  combination  to  derive  the 
calibration  coefficients  applied  to  all  flight  lines  on  a  single  day.  This  approach  does 
not,  however,  account  for  bidirectional  reflectance  differences  between  the  flight  lines 
(east-to-west  flight  orientation  vs.  west-to-east  orientation;  Figure  1),  so  some  banding 
artifacts  remain  in  the  final  products  at  the  junctures  of  the  flight  lines  (Figure  2, 
chlorophyll  mosaic). 

Prior  to  an  ELC,  we  first  selected  a  coastal  flight  line  with  a  ship  station  that  was 
collected  close  in  time.  For  example,  on  21  May,  we  extracted  the  spectral  PHILLS 
radiance  values  from  runl,  line  8  at  the  location  of  R/V  Ocean  Color  station  1  (see  Figure 
1).  This  flight  line  of  data  was  collected  from  0850  -  0856  local  time  and  the  ship  station 
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was  occupied  at  0858.  We  then  applied  the  ELC  technique  to  derive  calibration 
coefficients  to  convert  the  image  radiances  to  the  ship-measured  Rrs  values.  This  forces 
the  PHILLS  and  ship  Rrs  spectra  to  match  exactly  at  that  single  point  in  the  image.  These 
coefficients  were  subsequently  applied  to  each  pixel  in  the  flight  line,  and  to  each  flight 
line  for  that  day.  Comparisons  between  ELC-calculated  and  ship-measured  Rrs  spectra 
agreed  well  at  independent  stations  not  used  to  derive  the  coefficients  (Figure  3),  so  we 
adopted  the  ELC  approach  for  all  PHILLS  imagery. 


Figure  3.  Comparison  of  ship-measured  Rrs  (red  line)  and  ELC-derived 
PHILLS  Rrs  (green  line),  for  R/V  Ocean  Color  station  3  on  20  May,  2002. 


We  also  explored  the  effect  of  different  flight  line/station  combination  pairs  on  the 
calculated  calibration  coefficients  and  resulting  derived  optical  products.  On  three  of  the 
days,  several  flight  line/station  combinations  were  selected  for  ELC  to  develop  different 
sets  of  calibration  coefficients.  The  coefficients  were  subsequently  applied  to  the 
imagery,  and  the  agreement  between  PHILLS-derived  and  ship-measured  absorption  and 
scattering  coefficients  was  used  to  evaluate  the  coefficients  and  select  the  best  set.  For 
example,  on  21  May,  we  evaluated  two  other  potential  flight  line/station  combinations 
(run  5,  line  9  and  R/V  Ocean  Color  station  4;  run  9,  line  6  and  R/V  Ocean  Color  station 
5).  However,  the  errors  between  the  derived  and  measured  spectral  optical  properties 
were  larger  than  with  the  original  calibration  coefficients.  In  general,  the  different  flight 
line/station  combinations  used  in  the  ELC  did  not  have  as  large  an  impact  on  the  resulting 
optical  estimates  as  might  have  been  expected  (results  not  shown). 

NIR  SLOPE  ALGORITHM  (CALCULATING  a  AND  bb) 

Following  application  of  the  ELC  to  the  PHILLS  imagery,  a  NIR  slope  algorithm  was 
applied  to  the  Rrs  spectra  at  each  image  pixel,  to  derive  the  optical  properties  of  the  water. 
The  NIR  algorithm  is  based  on  the  observation  that  the  total  absorption  in  the  715  -  735 
nm  wavelength  range  is  controlled  by  pure-water  absorption  (i.e.,  the  absorption  by  the 
other  components  such  as  phytoplankton,  detritus,  and  CDOM  is  minimal).  Furthermore, 
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we  assume  that  the  spectral  shapes  of  b  and  bb  are  linear  and  relatively  flat  over  this 
narrow  wavelength  range  (only  a  2.8%  difference  in  b(715)  and  b(735),  based  on  the 
spectral  model  of  Gould  et  al.,  1999),  and  that  the  C  term  below  is  constant  (C  =  (t2*f)  /  ( 
n2*Q)  =  0.051).  These  assumptions  allow  calculation  of  the  backscattering  coefficient  at 
715  nm  from  the  corrected  Rrs  values,  and  subsequently  spectral  backscattering  and 
absorption  over  the  entire  400  -  700  nm  range.  The  calculations  are  summarized  below. 

Based  on  Rrs  =  C  bb  /  (a+bb)  (Gordon  et  al.,  1988)) 

Subscript  1  denotes  A.  =  715  nm,  subscript  2  denotes  A.  =  735  nm 
awi,  pure-water  absorption  at  715  nm  =  1.007 
aW2,  pure-water  absorption  at  735  nm  =  2.39 
bb2  =  0.97234  bb, 


First,  calculate  the  reflectance  difference  at  two  NIR  wavelengths: 

Rr.-<\  ~  Rrs2  ~  Cbbl  l{aw{  +  bbl)  —  Cbb2  /(aw2  +  bb2)  (1) 

Then,  substitute  and  rearrange  terms: 

0.972346M2  +  [0.97234^,  +  a„2  -C{aw2  -  0.97234a  J/^  +awlaw2  =  0  (2) 

Solve  the  quadratic  for  the  backscattering  coefficient  at  715  nm,  bbi: 

bbl  =  (-  B  -  ~Jb2  -9.362)/ 1 .945  (3) 

where  B  is  the  term  in  brackets  in  equation  2. 

Next,  calculate  spectral  bb(A,): 


bb  W  =  bbl  x  bbshape(A.)  (4) 

where  bbshape(X)  =  (1.62517  -0.001 13A)/ 0.8 17220  (5) 

from  the  spectral  model  of  Gould  et  al.  (1999). 

Convert  bb(A.)  to  b(A): 


b{X)=  53.56857  xbb  (A) +0.00765  (6) 

also  from  Gould  et  al.  (1999). 

Finally,  calculate  total  absorption,  aj(A): 

aT  (T)  =  Cbh  (A)/  -bb{X)  (7) 
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OPTICAL  WATER  MASS  CLASSIFICATION  ALGORITHMS  (PARTITIONING  a 
AND  b  INTO  COMPONENTS) 

Following  calculation  of  a,  bb,  and  b,  the  scattering  coefficient  is  partitioned  into  organic 
and  inorganic  components  and  the  total  absorption  coefficient  is  partitioned  into 
phytoplankton,  detrital,  and  CDOM  components,  using  empirical  relationships.  In 
addition,  we  require  estimates  of  chlorophyll  (here  we  use  the  Stumpf  et  al.  (2000) 
algorithm)  as  well  as  a<p(443)  derived  from  the  chlorophyll  value  (Stumpf  et  al.,  2003). 
The  algorithm  steps  can  be  summarized  as  follows: 

1.  Calculate  the  organic  component  of  scattering  at  660  nm,  bo(660),  from 
chlorophyll  (Loisel  &  Morel,  1998) 

2.  Calculate  b0(555)  from  bo(660) 

3.  Calculate  the  inorganic  component  of  scattering  at  555  nm,  bj(555)  =  b(555)  - 
b0(555) 

4.  Calculate  ad(4 12)  from  b;(55 5) 

5.  Calculate  total  absorption  minus  water  absorption,  ax-w(^) =  ax(A-)  -  aw(^) 

6.  Calculate  spectral  phytoplankton  absorption,  a^(X),  from  a,?(443) 

7.  Calculate  combined  detrital  and  CDOM  absorption,  a<jg(A,)  =  ax-w(^)  -  a,p(X) 

8 .  Calculate  CDOM  absorption  at  4 1 2  nm,  ag(4 1 2)  =  adg(4 1 2)  -  ad(4 1 2) 

9.  Calculate  %(A.)  from  ag(4 12) 

10.  Calculate  a<i(7.)  =  a<]g(A,)  -  ag(X) 

1 1 .  Calculate  particulate  absorption  ap(A,)  =  a,P(A.)  +  aj(A,) 

12.  Calculate  PIM  from  ap(4 1 2) 

13.  Calculate  POM  from  ap(443) 

14.  Calculate  TSS  =  PIM  +  POM 

See  Gould  et  al.  (2002)  for  equations,  references,  and  further  details  on  calculating 
organic  scattering  (b0),  ad(412),  ag(k),  PIM,  and  POM. 

RESULTS 

After  applying  the  ELC,  chlorophyll,  and  NIR  bio-optical  algorithms,  we  combine  the 
daily  flight  lines  to  create  image  mosaics  of  the  region,  for  each  of  the  optical  products. 
PHILLS  and  SeaWiFS  chlorophyll  images  covering  roughly  the  same  area  of  Mississippi 
Sound  near  Mobile  Bay  for  May  20  are  shown  in  Figure  4.  The  same  chlorophyll 
algorithm  was  applied  to  both  images,  and  the  color  scales  are  identical.  Although  the 
PHILLS  imagery  is  at  much  higher  spatial  resolution,  the  close  agreement  between  the 
two  images,  in  both  magnitude  and  pattern,  indicates  that  the  Rrs  values  resulting  from  the 
ELC  must  also  be  in  close  agreement  with  the  SeaWiFS  Rrs  values. 
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Figure  4.  Chlorophyll  concentration,  20  May,  2002.  A.  SeaWiFS.  B.  PHILLS 


Mosaics  of  organic  and  inorganic  scattering  coefficients,  CDOM,  detrital,  and 
phytoplankton  absorption  coefficients,  PIM,  and  POM,  for  21  May  are  shown  in  Figure 
5.  Note  the  different  scales  and  units.  The  scattering  coefficients  due  to  inorganic  matter 
are  in  general  an  order  of  magnitude  larger  than  the  organic  ones,  although  concentrations 
of  particulate  inorganic  matter  are  only  2-3  times  larger  than  concentrations  of  particulate 
organic  matter.  Also,  the  phytoplankton  absorption  coefficients  are  about  6  times  lower 
than  those  due  to  CDOM  and  detritus,  and  the  distribution  is  more  dispersed. 

Comparisons  of  PHILLS-derived  and  ac9-measured  spectral  absorption  and  scattering 
coefficients  from  R/V  Ocean  Color  station  1  on  21  May  are  shown  in  Figure  6.  This 
station  was  used  to  perform  the  ELC,  so  the  PHILLS  Rrs  spectrum  matches  the  ship 
measurements  exactly.  The  strong  agreement  between  the  PHILLS  and  ac9  values  at  this 
station  where  the  atmospheric  correction  was  “perfect”,  particularly  for  the  scattering 
coefficient,  indicates  the  efficacy  of  the  NIR  slope  algorithms.  A  non-constant  spectral 
shape  for  the  C  term  in  equation  1,  or  variability  in  the  bt/b  ratio  could  contribute  to  the 
discrepancy  between  the  derived  and  measured  absorption  values.  Calculations  of  error 
statistics  for  all  coincident  ship  and  aircraft  values  are  in  progress. 

We  use  the  partitioned  absorption  components  in  an  optical  water  mass  classification 
scheme  to  characterize  the  region.  The  percentages  of  total  absorption  due  to  ad(412), 
ag(412),  and  a<p(443)  form  the  three  axes  in  a  ternary  diagram,  to  quantitatively  illustrate 
the  dominant  absorbing  species  at  each  image  pixel.  Because  we  are  using  absorption 
percentages,  we  are  looking  at  relative  changes  in  the  dominant  absorbing  characteristics 
at  each  pixel  rather  than  changes  in  absolute  magnitude  of  the  absorption  coefficients. 
Thus  the  patterns  show  compositional  changes.  The  optical  classification  scheme  allows 
us  to  examine  spatial  and  temporal  optical  variability,  as  well  as  follow  features  and  track 
water  masses.  Figure  7  shows  an  example  of  this  approach  applied  to  the  19  May 
PHILLS  imagery.  Red  pixels  in  the  RGB  composite  mage  indicate  relatively 
higherdetrital  absorption,  green  pixels  indicate  relatively  higher  phytoplankton 
absorption,  and  blue  pixels  indicate  relatively  higher  CDOM  absorption. 
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Figure  6.  Total  absorption  and  scattering  coefficients  vs.  wavelength, 
PHILLS-derived  and  ac9-measured,  R/V  Ocean  Color  station  1  on  21  May. 


Figure  7.  Optical  water  mass  classification  for  PHILLS  imagery,  19  May.  RGB 
composite  image  with  %ad(412))  loaded  in  the  red  channel,  %a<p(443)  loaded  in  the 
green  channel,  and  %ag(412)  loaded  in  the  blue  channel.  The  black  points  clustered  in 
the  lower  right  corner  of  the  ternary  diagram  indicate  that  although  there  is  some 
variability,  CDOM  absorption  dominates  the  total  absorption  signal  in  this  region. 
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DISCUSSION  AND  SUMMARY 


We  have  presented  new  algorithms  to  estimate  backscattering  and  absorption  coefficients 
from  hyperspectral  ocean  color  imagery;  based  on  spectral  RrS  at  NIR  wavelengths. 
Empirical  relationships  are  used  to  partition  the  absorption  coefficient  into 
phytoplankton,  detrital,  and  CDOM  components,  and  these  components  are  ultimately 
used  in  an  optical  water  mass  classification  scheme.  These  approaches  were  applied  to 
high-resolution  (10  m  pixels),  hyperspectral  aircraft  imagery  collected  near  Mobile  Bay, 
Alabama.  Several  sources  of  error  contribute  to  the  overall  agreement  between  PHILLS 
estimates  and  ship  measurements: 

1.  PHILLS  sensor  calibration 

2.  Empirical  Line  Calibration  method  -  incorrect  RrS  estimates 

3.  Non-negligible  CDOM,  detrital,  or  pigment  absorption  at  NIR  wavelengths 

4.  C  term  incorrect  or  not  constant 

5.  Incorrect  shape  of  spectral  b(X)  or  bt,(X) 

6.  Incorrect  or  variable  bb(A)/b(Z)  ratio  (i.e.,  variable  volume-scattering  function 
between  stations) 

We  are  exploring  these  error  sources  further,  but  initial  results  show  good  agreement 
between  PHILLS-derived  values  and  ship  measurements  (both  Rrs  and  optical  properties). 
The  optical  water  mass  classification  technique  (RGB  imagery  and  ternary  plots)  helped 
delineate  patterns  in  the  absorbing  characteristics  of  the  waters  in  and  around  Mississippi 
Sound,  and  indicated  that  CDOM  dominates  absorption  in  this  region. 
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